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Abstract. The evolution of three-dimensional disturbances in an incompressible three- 
dimensional stagnation-point flow in an inviscid fluid is investigated. Since it is not possible to 
apply classical normal mode analysis to the disturbance equations for the fully three-dimensional 
stagnation-point flow to obtain solutions, an initial-value problem is solved instead. The evolution 
of the disturbances provide the necessary information to determine stability and indeed the com- 
plete transient as well. It is found that when considering the disturbance energy, the planar 
stagnation-point flow, which is independent of one of the transverse coordinates, represents a neu- 
trally stable flow whereas the fully three-dimensional flow is either stable or unstable, depending 
on whether the flow is away from or towards the stagnation point in the transverse direction that is 
neglected in the planar stagnation point. 
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1. Introduction. 


Previous analytical work investigating the stability of planar stagnation-point flows has con- 
centrated on the mathematical simplification provided by classical mode analysis of streamwise dis- 
turbances. Wilson and Gladwell (1978) have shown that incompressible planar stagnation flow is 
always stable to three-dimensional normal mode self-similar disturbances that decay exponentially 
outside of the viscous boundary layer. Lyell and Huerre (1985) re-examined the planar stagnation 
flow problem by using the same class of disturbances and verified the results of Wilson and 
Gladwell as well as characterized the other stable eigenvalue branches by showing that, after the 
initial branch found by Wilson and Gladwell, the other branches come in pairs. In addition, by a 
nonlinear analysis using a Galerkin method, these authors indicated that this same flow is unstable 
for disturbances of sufficiently high amplitude. On the other hand, a numerical study by Spalart 
(1989) found no such instability suggesting that this was not the case. Brattkus and Davis (1991) 
showed that the normal mode self-similar disturbances were the least stable of the class of distur- 
bances that have a power like behaviour in the downstream coordinate. 

Lasseigne and Jackson (1992) allowed for density variations induced by a temperature 
difference between the freestream and the plate and determined that the stagnation flow remained 
stable to small streamwise disturbances regardless of the plate temperature. The effect of cooling 
the plate was to decrease the decay rate (less stable) of the small wavelength disturbances while 
increasing the decay rate (more stable) of the moderate wavelengths. Again, only three-dimensional 
normal mode self-similar disturbances that decay exponentially outside of the viscous boundary 
layer were considered. Studies dealing with the swept attachment line also concentrated on stream- 
wise normal mode linear disturbances. Hall, Malik and Poll (1984) determined that a region of 
instability (in frequency-wave number space) associated with increasing crossflow exists. In 
independent investigations Kazakov (1991) and Lasseigne, Jackson and Hu (1992) determined the 
effects of surface temperature variations on this region of instability; the latter investigation also 
allowed for the effects of suction or blowing at the surface. 

It is not possible to analyze the stability of the fully three-dimensional stagnation-point flow 
using normal mode streamwise disturbances since the disturbance equations do not admit this class 
of disturbances as an eigenvalue problem. Therefore, the approach taken in this investigation is 
different than that of previous investigations. A more general initial value problem is solved using 
the methodology developed in Criminale and Drazin (1990) and has its origins from the work of 
Kelvin (1887) and Orr (1907a,b). The disturbances are taken to be initially bounded in all direc- 
tions and the evolution of initial conditions is determined analytically and in closed form. By con- 
centrating on the mean flow subject to disturbances in an inviscid fluid, the fully three-dimensional 
stagnation-point flow can be solved with the planar stagnation-point flow as a special case. 

The method of analysis utilizes a moving coordinate transformation that allows for easy 
integration (in time) of the individual vorticity components. Then, the double Fourier transform in 
the new transverse coordinates is used to reduce the mathematical problem to the solution of ordi- 
nary differential equations in which time appears strictly as a parameter. Thus, a completely analyt- 
ical solution is found to the initial value problem describing the evolution of three-dimensional dis- 
turbances in three-dimensional stagnation-point flow in an inviscid fluid. The time evolution of a 
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single Fourier mode in which the disturbance is periodic in the transverse directions is investigated 
in detail as well as the evolution of an initially localized disturbance. In both cases, the evolution 
of the total energy of the disturbance is used to illustrate the importance of the transient. Kelvin 
(1887) and Orr (I907a,b) have shown for different flow configurations that there can be growth in 
the perturbation energy even when the classical mode analysis shows that the flow is stable for 
long times. 

Farrell (1989) applied this approach to two-dimensional disturbances within a two- 
dimensional counterflow as a good approximation to the local flow in regions of confluence and 
diffluence. It was determined that plane wave disturbances with dependence in the transverse coor- 
dinate and independent of the other coordinate have energy that grows exponentially in time. Plane 
wave disturbances with dependence in both coordinates experience an initial exponential growth, 
but eventually decay as time progresses. Disturbances with finite wave trains that were not spatially 
symmetric were shown to be stable but these disturbances also experienced an initial transient 
growth in energy with the energy asymptotically approaching a constant amplitude. Symmetric 
finite wave train disturbances were shown to not experience the initial transient growth and to have 
energy constant in time. This was seen to occur in both regions of confluence and diffluence. 

The governing equations for the three-dimensional stagnation-point flow, the moving coordi- 
nate transformation and the linear disturbance equations are presented with the disturbance equa- 
tions solved by the use of Fourier transforms in the transverse coordinates. Selected results for the 
evolution of a single Fourier mode are given in Section 3 and, in Section 4, results pertaining to 
the time evolution of a finite wave train described by an initially Gaussian profile are presented. In 
Section 5, a constant pressure boundary condition is considered as an alternative to the zero normal 
velocity condition. Section 6 contains a discussion on the effects of background rotation and parti- 
cle paths. Conclusions are given in Section 7. 

2. Problem Statement and Basic Equations. 

The problem under investigation is that of linearized disturbances in a three-dimensional 
stagnation-point flow. The basic flow is given nondimensionally by 

U = x, V =-(l+X)y, W = Xz, (1) 

where X is a measure of three-dimensionality. For X = 0 the flow is a two-dimensional stagnation- 
point flow; for X = 1 the flow is axisymmetric; for - 0.4294 < A, < 0 the flow corresponds to two 
symmetrically displaced protuberances (Davey, 1961) displayed in Figure I. For X < 0, it is 
important to observe that the flow is toward the stagnation point in the z -direction and away from 
the stagnation point in the x -direction. In addition, Davey has shown that separation occurs at 
X = -0.4294, reversed flow exists for - 1 < X < -0.4294, and no solutions are possible for 
X<-\. 

The nondimensional linearized equations for small disturbances are written as 

U X + Vy + W Z = 0 


(2a) 



(2b) 
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u ( + x u x - (l +\)y u y + \z u z + u = -p x 


v, + x v x - (1 + X)y v y + Xz v z - (1 +X) v = -p y 


(2c) 


w, + xw x - (1 + X)y w y + Xz w z + = -p z 


(2d) 


where w,v,w and p are the velocity and pressure perturbations, respectively. The appropriate 
boundary conditions require that all disturbance quantities vanish as y -> °° and the normal velo- 
city v is zero on the wall. In addition, initial conditions consistent with the boundary conditions 
must be supplied. 


The above equations can be recast in terms of the vorticity components ciy, , co v , w z in the 
x,y,z directions and are 


Dt 


= (0, 


D(O y 

Dt 


-(1 + X)co ) , , 


D (O z 
Dt 


= X(l) z , 


(3) 


where D IDt is the material derivative defined by 

= ( )< + x ( )j( — (1 + X,) y ( ) y + X,z() z . 


(4) 


In general, only two vorticity components can be specified at time t = 0; the third component is 


found by appealing to the equation of conservation of vorticity, given by 

9(o, 9(o v 3<n z 

-af + Tf + l r= 0 ' 


(5) 


In this study, the initial profiles for the vorticity components <0* and (o z are specified, and the 
above equation is used to determine (0^. Once the vorticity has been determined from (3), the v 
velocity component is found from the following relationship 


V 2 v 


9ca z 

dx 


9% 

9z 


( 6 ) 


where V 2 is the three-dimensional Laplace operator. The other velocity components are determined 
by appealing to the vorticity relations together with the continuity equation, yielding: 



9 2 v 

9x9y 


(7) 


V 2 w = - 


9co y 

rh 


9 2 v 

9y9z 


( 8 ) 


where is the two-dimensional Laplace operator in the x - z plane. 
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By following Criminale and Drazin (1990), an analytical solution is readily obtained if a 
change from an Eulerian to a moving coordinate system is made. The moving coordinate transfor- 
mation 

^ = xe~', r|=yc (1+X) ', C, = ze~ Xl , T = t (9) 


is chosen so that the system of equations (3) have coefficients that are functions of time only and 
the material derivative (4) becomes 


DO _ 90 
Dt dT ' 


( 10 ) 


The notable feature of this transformation is that the partial differential equations for the individual 
vorticity components (3) can then be immediately integrated in time and the solutions are 


0 T 

(0 X = <tt,V 


CO, = co,V (,+X)r 


®z = ® z V r 


(ID 


where the superscript 0 denotes the specified values at T = 0. Once the initial values co° and o) z ° 
have been made, the initial value of co° is determined from the conservation of vorticity (equation 
(5)). The equation for the v velocity component in the moving coordinates is 


Av = e~ (I ~ X)T 3<0z ° - e°~ X)T d( °* 


* 




( 12 ) 


where A is the Laplace operator in the moving coordinates, namely 


A = e 


-2 T 


as 2 


+ e 


2(1 +A.)T 


an 2 


+ e 


-2 XT 


ac 2 


(13) 


It is noted that the time dependence in equation (12) appears as a parameter only and hence finding 
a solution to (12) is essentially a spatial problem. The equations for the other two velocity com- 
ponents in the moving coordinate system are 


A 2 « 


I).,. 0 

_ -M + 2X1T pu 7 _ \T 

ac 


a 2 v 

acan 


(14) 


A 2 w 


-<2-i-x)r a®> _ t a 2 v 

ac acan ’ 


(15) 


where 



+ e 


-2 XT 



(16) 


is the two-dimensional Laplace operator in the moving system. 

It is now assumed that Fourier transforms may be taken in the C and C directions, which 
implies that the disturbances are bounded in these directions. The double Fourier transform is 
defined by 


& 


\ 



oo 


(17) 
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«(a,r|,Y, T) = JJ u(Z,,r\,{ 3 ,T)e l(a t> + yl ' ) dt,d( ) , 

— oo 


etc. for <?, w, p, d) with the inversion given by 

oo 

u & ri, C, T) = JJ u (a, ri, y, T) e' 1 <°* + ™ dady. 

(lit) -eo 


(18) 


The transforms of the relevant equations are 


36° 

5r| 


= i at})® + i yti>z , 


(19) 


( 20 ) 



i y&y e 


0 < ,-(l+2X)T _ ia jW_ p >-T 




( 21 ) 



.fad>,V‘ 2+w -iy^-e T 
y 5t1 


with 

a 2 = a 2 e~ 2T + y l e~ 2XT . 


( 22 ) 


(23) 


The equation for v is a second order ordinary differential equation in rj with time dependent 
coefficients, and since ft 2 is just the Fourier transform of A 2 , the velocity components u and h> are 
determined directly from the algebraic relations provided by transforming equations (14) and (15). 

The evolution of a single Fourier mode can be studied by choosing the initial conditions for 
the vorticity components to* and co z ° to be 

(w°,co z °) = (Q* ,£2 Z )e~'' (a °^ + To?) 5(Ti-y 0 ) (24) 

since (x,y,z) = Q at T = 0. In Fourier space, 

(<b°,fi> z °) = (Q, ,n z )5(a-oto)8(y-Yo)5(ii-yo)- (25) 

The initial value for the last vorticity component is found by transforming equation (19) and 
integrating in the rpdirection: 

d) y ° = (/ an x + i y£2 z ) 5(a - ao) 5 ( y - y 0 ) H(r\ - y 0 ) - 1 j 


(26) 
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where H (r|) is the step function. The constant of integration has been chosen so that all of the 
velocity components vanish as q — > 

With the vorticity components known, the differential equation for v is integrated, and 
V = - ^ ^ — | / aL2 z e -^ + ^ T - j yQ x | §^ a _ a 0 )5(y - y 0 ) 

x e -tt^' 0+x,r l n-yol _ g -ft‘’' (1+x,T (n+yo)j ^27) 


where Ot is as defined earlier. The solution decays as q — » °° and is zero at the wall. Once the 
transform is inverted, the velocity component is given by 

= tAt j / - ,y„ll x 

2(2 Jtr&o [ J 

x ,-<V-‘ ,U,r N->ol _ -*o‘-«"> r < n +y4] ^28) 


where 


a 0 2 = a, 2 e“ 21 ' + ySe~ 2lT 


The u and w components arc determined directly from the earlier equations and, upon inverting 
the transforms these components are found to be 


H&r|,C,r) = 


( 27 t ) 2 &0 


2(2ti) 2 & 0 2 


<*04, + 7 o 4 ] ^ (,+2X)r e' _, ' (ao ^ +Y ° g [ H( n - y 0 ) - 1 


otoQ, e~ JT -y 0 n x e-^ T i 


x [ sign (t 1 -y o) r ftoe ' < " X,r| ^ 


w(^,%^T)= \ Oofl, +Tfo4rl e- (2+X)7 V' (a °^ Yo0 [ H (r\ — y$) — 1 

(2tt) 2 &o L J L 


-\-j \ a«n z e-^ T -y 0 n x e~^) e- 1 ^^ 
2 ( 27 t) 2 &o 2 [ J 

x ’ sign m-y 0 ) e -^ ( ' +X)Tl ^ - e -*o-- a "> T *+y*] . 


3. Single Mode: Infinite Wave Train. 

In this section we present selected results for the single mode solutions given above for u, v , 
and w . We first note that the solutions are a linear combination of Q x and 4 • Therefore, there 
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exists two modes, defined as follows: 

Mode / : £2, = 0, £2 Z = 1 

Mode II : £\ = 1, £2 Z = 0 


The results for Mode I are shown in Figure 2 for three values of the wavenumber in the z - 
direction (y 0 = 0, y 0 = 0.1, and y 0 = 1.0) and a single value of the wavenumber in the ;t -direction 
(ocq = 1.0). Plotted in this figure is the maximum amplitude of the velocity components I u I maz , 

I v I max’ and I w I max as a function of time. The maximum in the u and w velocities are at the 
wall, while the maximum in the v velocity is in the interior for all times since it is required that v 
vanish at the wall. In each sub-figure, four values of the three-dimensionality constant were 
chosen: a) axisymmetric flow, X = 1 (shown in the figures as a small-dashed curve), b) X = 0.5, c) 
planar stagnation-point flow, X = 0 (shown in the figures as a large-dashed curve), and d) 
X = -0.2. 

Included in these figures are the results for the special cases of two-dimensional disturbances 
in a two-dimensional flow (yo = 0, X = 0), two-dimensional disturbances in a three-dimensional 
flow (y 0 = 0, X * 0), and three-dimensional disturbances in a two-dimensional flow (y 0 * 0, X = 0). 

As seen from the figures, all three components of the flow exhibit exponential decay in the 
long time solution for the response of this mode. The results for the u velocity show that the decay 
rate is independent of die variable X and of the wave number in the z direction. Since the analytic 
solution for the velocities is known, this information can be found by determining the behaviour of 
equation (30) as T — » « and is 

M(t0,C’F)= — — e - r £'‘ i(ao5+Yo0 (32) 

(2 it) 2 

where the first term of (30) is dominant when X < 1, and there is a balance between the first and 
second terms for the axisymmetric case X = 1 that does not, however, change the limit. For the 
long time behaviour of the v component of the velocity the maximum must be determined. The 
maximum is found to always occurs at q = y 0 so that 

v&yo.C.n = e -(W)r (33) 

(2tc) 2 

with the axisymmetric case again being a special limit that does affect the amplitude in the above 
equation but not the decay rate. The asymptotic behaviour for the w component of the velocity is 
also determined analytically. From equation (31) it is found that 

w&O ,;.T) (34) 

(2 It) 2 I Yo I 

The dependence of the decay rate of the v and w components on the three-dimensionality parame- 
ter X is clearly seen in Figure 2. 


- 8 - 


Although die above analytic forms appear to be relatively simple, the use of the transformed 
variables t, and r| tend to hide some rather significant changes in the spatial structure of the dis- 
turbances. For X, > 0 there is a stretching of the initial Fourier mode form in both the x -direction 
and the z -direction. For the axisymmetric case, this stretching occurs at an equal rate in both direc- 
tions. For the planar stagnation-point flow X = 0 there is no stretching of the disturbance in the z 
direction, and for X < 0 where the mean flow is towards the stagnation point in the z -direction, 
there is an actual contraction of the initial Fourier mode form in the z -direction. In the y direction, 
it is noted that although the maximum in v velocity is at a fixed value of r), this implies that this 
maximum approaches the wall exponentially fast when the problem is converted to the physical 
spatial variables, that is the inverse of (9). 

Although the long time behaviour of the disturbances is of course very important, equally 
important is the transient evolution since it is possible for there to be significant initial growth (or 
perhaps a long time persistence) of the disturbance quantities before the inevitable exponential 
decay. By examining the results of Figure 2, it is seen that the initial response to a Mode I distur- 
bance is typically a linear decay in time. The initial linear decay rate depends on the value of the 
three-dimensionality parameter X with X < 0 decaying at the slowest rate. The X dependence of this 
initial decay rate is strongest for smaller wave numbers in the z -direction. 

The results for the single mode response to a Mode II disturbance are shown in Figure 3. A 
significant difference from the Mode I results is immediately noticeable. Analytically the long-time 
limits are given by 

„(l.o ,;,n= “°T° >, ‘ 1 , (35) 

(2 k ) 2 I Y 0 1 


v&yo .c.n = - 


[ Yoyo 
(2 k ) 2 


-(l+3X.)r (“oS+YoO 

t c j 


(36) 


and 


w&o,c,r) = - 


i 

(2 k ) 2 


\T ^-'(OoS+YoS) 


(37) 


as 7 1 ->«, it is seen from Figure 3 that the m and v components of velocity decay exponentially 
for all X while the w component of velocity only decays when the three-dimensionality parameter 
X is greater than zero. This value of X corresponds to a stretching of the original disturbance in the 
z direction. However, for the planar stagnation-point flow where X = 0 the w component 
approaches a constant independent of the initial wave number in the x and z directions. For the 
case where there is mean flow towards the stagnation point along the z -direction (X < 0), the w 
component grows exponentially. It is seen that the initial behaviour of the u and v components 
also reflect this differing behaviour for the planar stagnation-point flow and the X < 0 case. It is 
also true that there is initially linear growth of these components. As in the Mode I case, the 
dependence of the behaviour for early time on the three-dimensionality parameter is strongest at 
smaller wave numbers in the z -direction 
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The results presented in Figures 2 and 3 give the maximum values of the velocities, but con- 
sidering the severe contraction of the disturbances in the y -direction as evidenced by the movement 
of the maximum in the v velocity toward the wall at an exponentially fast rate, it is perhaps more 
useful to look at the growth or decay of the energy of the disturbances. For the single Fourier 
mode results presented in this section, it is necessary- to consider the energy per period in the \ and 
£ directions, or equivalently the energy per period in the x and 2 directions. Because of the 
stretching in the x -direction and the stretching or contraction (depending on the sign of X) in the 
2 -direction, these results do not readily lend themselves to direct physical interpretation; however, 
they do help point out the effect that the contraction has in the vertical direction y. Shown in Fig- 
ure 4 is the quantity £ as a function of time defined by 




l 2 + I v I 2 + Iwl 2 


dy 


(38) 


where w, v, and vc are given by (28), (30), and (31). Each graph was normalized by the value at 
time / = 0. The amplitudes and Q 2 are equal in the mixed mode results. In each sub-figure, 
four values of the three-dimensionality constant were chosen: a) axisymmetric flow, X = 1 (shown 
in the figures as a small-dashed curve), b) planar stagnation-point flow, X — 0, c) X — — 1/3 
(shown -in the figures as a large-dashed curve), and d) X = -0.4. From these figures, Mode II is 
seen to, represent disturbances that may grow in energy as it also represented disturbances with 
growth in the maximum values of the velocities. However, it is noted that although the maximum 
in the velocities can grow for any X < 0, the energy per period grows only for X < — 1 / 3. This 
phenomenon is directly related to the contraction in the y direction. Very little dependence on the 
wave number in the 2 direction can be detected which is consistent with the results of Figures 2 
and 3. The energy of Mode 1 is dominated by the u and v velocities which show little depen- 
dence on y 0 , and the Mode II energy is dominated by the w velocity which also shows little depen- 
dence on y 0 . 


4. Finite Wave Packet. 


In the previous section, the results for an initial disturbance that consisted of a single Fourier 
mode were presented. This disturbance has an infinite spatial profile in the a - 2 plane. Considering 
the distortions in the x and 2 directions introduced by the transformed variables, there is consider- 
able doubt as to the proper interpretation of the single mode results. To resolve this difficulty in 
the analysis, an initial disturbance that is localized in space is chosen. The initial vorticity profile is 
taken as a Gaussian. This disturbance is initially symmetric in the x -2 plane and, since the invis- 
cid problem is being investigated, a Gaussian profile in the y direction is also chosen and all boun- 
dary conditions can still be satisfied. 


The analysis proceeds by replacing (24) with 


((0°C0.°) = 


(Q,.flz) 

(2 7 tV [}) ,/4 


_ii +£ __al 

4 a 4 B 


e 


(39) 
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In Fourier space, 


«&°.<b- 0 ) = V °(a^rV 4 

(2 7 r cr P ) /4 


_nl 


(40) 


where the constant is chosen such that the integral of the square of the vorticity over the domain is 
a constant for each choice of o and p. 

The above initial conditions can be substituted directly into equations (19) and (20) which 
can then be solved for l? and <5^° or the previous results for a delta function initial condition in the 
y direction can be integrated over the dummy variable y 0 . Either way, the transform of the vorti- 
city (0° and the v component of velocity are given by 


<£>,° = -(' + i yR) 4 erfc 

(2 J ro‘ P) 




(41) 


and 


v =(/ an z c- (?+l)T - i yQ x e~ {urK)T ) 4n - 

z x (2tt , o 2 P) ' 4 


X 1 e $Ci 7 c- 2 < uX)T e (i+X)T 
2 & 




erf c 


erfc 


_JL_ &t> -o + >orvp 

2 Vp H 


+ 2e 


Uc 


_2L + &t ,-(t + X)rvp 

2Vp 


°’ x,r Vr/ [&e- (1+X)r V p 


■ (42) 


The transforms of the u and w< components are found directly from equations (21) and (22) and all 
transforms can be inverted using the formula (18). The long time behaviour of w at the wall for 
Mode II is given by 


(0, T) = - — 47iq . ffi* _ e -<na^f) e -^T 

(2k 3 o 2 P) ,m 


(43) 


as 7" — > oo which shows that the transform grows exponentially for X negative consistent with the 
results of the single mode. 

In order to examine the true behaviour of the evolution to this initially localized disturbance, 
the energy as a function of time is computed. The energy in terms of the unstretched variables is 




E (T) = — | II [w 2 + v 2 + w 2 ]dx dz dy , 


(44) 


which can also be given in terms of integrals over the Fourier transform quantities in the stretched 
variables as 


E(T) = 


,d+x)r 

Tn>~ 


Q\&e T ,ye XT ,y e (,+X)T ,T)dHdydy , 


(45) 
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where 

<2 2 (a,Y,rt,T) = « 2 +v 2 +vv 2 . (46) 

The normalized energy of the response to a Mode II gaussian profile with P = 1.0 and o = 0.5 is 
shown in Figure 5 for eight values of the three-dimensionality parameter X ranging from X = 1 (the 
axisymmetric case) to X = — .4 (near separation). A number of interesting features are noticeable. 
First, for the planar stagnation-point flow (X = 0), the energy approaches a constant value approxi- 
mately three times the initial value. This is in agreement with the results tracking the maximum in 
the velocities. The mean flow produces a coniraction of the disturbance in the y -direction, initially 
increasing the magnitude of the w velocity at the wall and initially increasing the energy of the 
disturbance. However, a balance develops between the continued contraction in the y -direction and 
the relief provided to the flow in the x -direction thereby leading to the situation where the energy 
approaches a constant. This behaviour of the disturbance energy for the two-dimensional planar 
counter-flow was also seen by Farrell (1989). The purpose of this study is to determine the effects 
of the three-dimensionality of the mean flow. When these three-dimensional effects are considered 
(X * 0), it is seen that the planar stagnation point flow represents a special case. For X slightly 
positive, there is again an initial transient growth in energy. The extra relief provided to the flow 7 
in the z -direction means that the balance that developed for the planar stagnation-point flow 7 does 
not develop and the energy of the disturbance eventually decays to zero. For X sufficiently large 
and positive there is no initial transient growth implying that the relief provided to the flow 7 in the 
z -direction is sufficient to prevent the the initial increase in disturbance energy. From Figure 5, it 
is seen that for the special case of axisymmetric flow (X = 1), the energy undergoes an immediate 
exponential decay. For X less than zero, an unstable situation develops, and the energy of the dis- 
turbance continually increases at an exponential rate owing to the contraction in both the y and z 
directions which is not balanced by the relief in the .r direction. 

It is also interesting to examine the behaviour of the vorticity components which are given by 
equation (11) in the transformed coordinates. In the physical coordinates, the vorticity field is 
undergoing a contraction in the y -direction, an expansion in the x -direction and an expansion (for 
X > 0) or a contraction (for X < 0) in the z -direction. The exponential time factors in equation (II) 
indicate that for a Mode II disturbance vorticity is transferred from tire y -component to the x- 
component while die transfer of vorticity between the y and z components for a Mode I distur- 
bance depends on the sign of X. By considering the response of a Mode II initial disturbance, it is 
seen that in the neutrally stable case (X = 0) the decay rate of the y -vorticity component is equal to 
the growth rate of the x -vorticity component. For the stable case (X > 0), the decay rate of the y- 
vorticity component is greater than the growlh rate of the a -vorticity component; whereas, for die 
unstable case (X < 0), the decay rate of the y -vorticity component is less than die growlh rate of 
die x -vorticity component. For the Mode I disturbance, either both components are decaying 
(X < 0), or the decay rate of the y -vorticity component is always greater than the growth rate of 
the z -vorticity component. Perhaps diis is an indication of why it is the Mode II disturbances that 
experience unstable energy growdi rather than those of Mode I. 
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5. Constant Pressure Boundary Condition. 

One problem of interest in geophysical fluid dynamics is the case where there is a constant 
pressure surface in lieu of a condition on the velocity. Such a problem has been examined by Eady 
(1949), for example, where plane Couette flow is mathematically equivalent when the constant 
pressure surface is used. Consider p = 0 at y = 0 rather than the v velocity. Of course the surface 
y = 0 can no longer be thought of as a solid wall. To derive an equation which governs the time 
evolution of the pressure, the momentum equations (2b) and (2c) can be combined with the con- 
tinuity equation (2a) in the (£,r|,£, T) coordinate system to get 

A 2 p =2a-l)<- r J| + r“ r ^<e lwi)T |^ (47) 


where it is immediately seen that the axisymmelric case X. = 1 is a special case. In keeping with 
the previous analysis, the pressure equation is first transformed into Fourier space and then u is 
eliminated by using (21), yielding 


.2 
& p = 


= -UsLa-i)&°e- 2 ^ l > r + a^i(x-i )e -('-»7 i£ 


or 




dr\ 


.,-2X7 d , (i + 3X)7 dv. 

dr 1 dn 


(48) 


If the initial conditions are considered to be given by equation (24) for this problem so that only a 
single Fourier mode is investigated, then is given by (26), while (19) can be integrated to find 
v, yielding 


v = A l (.T)e- hc ~ (UhT]r '~ y ° l +A 2 (T)e &c ‘ lU? ' ,r 


(49) 


where 


A X (T) = 


,U +X)T 
2 & 


i<xn.e~ l7+X)T -i yn x e~ (U7X)T [ 6(a - exo) 6(y - y 0 ) . 


(50) 


The determination of A ,(T) satisfies all of the conditions at r| = y 0 and it is left to impose the con- 
dition on the pressure at T| = 0 in order to determine the remaining unknown function A 2(7"). 
Equation (48) provides the necessary information in determining A 2 , but it is convenient to work 
with an alternate dependent variable defined by 


B(T) — e {Un)T -^(0,T) = he 2XT 
dr) 






•d -hr 


— A 2 V 


& ynC 


-d + 


(51) 


The equation for B(T) is then 


&- lsi ( x-i) P -77 g= _im ( x-i)n f - ir 

dT & 2 & 2 


(52) 


where 
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ft = (/aft, +i yt2 z )8(a- Oo)5(y- Yo). 


(53) 


It is helpful to find the long time solution in order to determine the effects of the three- 
dimensionality parameter X. Since Mode II (ft* = l,ft 2 =0) is the more interesting initial condi- 
tion, the rest of the analysis is restricted to this case with X< 1. The long time solution for B is 
determined directly from (52) or 

B =B„U- , (54) 

r 

as T —> oo where the the proportionality to ft is explicitly given and the coastant B^ can be found 
numerically. The behaviour of A ( is determined directly from equation (50), and A 2 from equation 
(51). These limits are 

A i(T) = -yr -r^r- e~ XT 5 (« " «o> S (Y - Yo) - (55) 

Z I Y i 

A 2 ( T ) = --TT^ -aBJ<?" xr 5(a-a 0 )5(Y-Yo)> (56) 

I y I 2 

as T — > oo. The long time behaviour of the velocities u,v,w is found through equations (24), 
(49) and (22) indicating 

u (0, T) = 4 (Y- «B») e~ T 8(a - cto) 5(y - y 0 ) , (57) 

T 

$ (0,T) = --r~r(y- aBJ e~ XT 5(a - 6(y - y 0 ) , (58) 

lYl 

w(0,T)=- — B„ e- XT 8(a - «o) 5(y - y 0 ) , (59) 

Y 

as T — > 00 . It is seen that, for the boundary condition p = 0 on the surface y = 0, the velocity 
components v and w of Mode II are exponentially growing when the three-dimensionality parame- 
ter X is negative. 

6. Effects of Background Rotation; Particle Paths. 

The expression used for the basic velocity and given by (I) can be modified to include effects 
due to background rotation (strain). Specifically, if ft 0 is the constant dimensional rate of rotation, 
then 

U= Jc+ft 0 z, y=-(l + \)y, W — Xz — ft„ jr , (60) 

becomes the new non-dimensional representation. The process of shifting to a moving coordinate 
system and solving for the perturbations as previously done can likewise be accomplished here but 
with a noticeable increase in complexity. For example, the fundamental vorticity equations (3) used 
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for solving the initial-value problem are no longer uncoupled and the net effect is a higher order 
differential system. Thus, the dynamics is altered to some degree. In particular, temporal oscilla- 
tions become possible but the overall stability conclusions remain. These points can be illustrated 
be examining the paths of the material particles in the basic flow. 

The Cartesian components for any particle in a Lagrangian frame can be written using (60) as 


— = u= x+ n fl z, 
dT 0 

(61a) 

iL = V=-(l +X)y, 
al 

(61b) 


- — = w = Xz - n p x . 
dT 


By assuming x = x 0 , y = y 0 and z — z 0 at time T = 0, the solutions for (61a,b,c) are 

x = x 0 e cT , 

y=y^- (1+x>r . 


(61c) 


(62a) 

(62b) 


1 


XT 


f p aT _ XT ' 

oX L J ’ 


(62c) 


with 


a = 


1 + X 
2 


1 + 


(1 + , 


1/2 "t 


(63) 


Oscillatory solutions are possible if > y (1 - X). In general, particles will move arbitrarily far 
from any initial position. 


Two interesting limits of (62a, b,c) are when n o = 0 or X = 1; X can never be very negative 
and therefore a change in the system behaviour using this parameter is not possible. In the first 
instance, = 0, then 


y =y a e 


-( 1 +X)T 


,XT 


(65) 


The novel feature of (65) is that these results are identical to those of (11) for the vorticity where 
x,y,z are replaced by to* , (o y , co z . In short, this is a unique situation where, for a three- 
dimensional flow, the particle paths are synonymous with vortex lines. Extending this argument to 
the n o * 0 case can be made by conjecture since the vorticity has not been determined under these 
circumstances. There is little likelihood that the unstable aspects of the flow, however, can be 
changed by a finite . When X = 1 , the basic flow is axisymmetric and o = I + / £2 C indicating 
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that background rotation has caused the system to have inertia and therefore oscillations in time. 

7. Conclusions. 

We have investigated the evolution of three-dimensional disturbances in a fully three- 
dimensional stagnation-point flow in an inviscid fluid. It has been shown that the planar 
stagnation-point flow is a special case in which the disturbance energy approaches a constant for 
long time. If the flow in the second transverse coordinate is away from the stagnation point then 
the flow is provided enough relief such that the disturbance energy decays after an initial transient 
growth. As the limiting case of axisymmetric flow is approached, the disturbance energy is found 
to decay without an initial transient growth. For flow towards the stagnation point in the second 
transverse direction, it is found that the disturbance energy may grow exponentially thus indicating 
an unstable flow configuration. These results were found by determining a closed form solution to 
the initial value problem even though a classical mode analysis was not possible for the fully 
three-dimensional flow. 

Because of the inviscid assumption, the results for the planar stagnation-point flow cannot be 
compared directly with previous work but the method used here can be extended (with consider- 
ably more mathematical complexity) to the study of the inviscid mean flow subject to viscous 
linear disturbances. These results can be compared with some of the previous work. However, any 
results for this problem can only be suggestive in view of the fact that the basic flow is derived 
from an inviscid basis. 

Background rotation or strain of the field can only alter the dynamics by forcing the system 
to have temporal oscillations rather than prevent instabilities. Finally, more comprehensive material 
particle path information could be obtained if one allows for the velocity field to include the pertur- 
bations as well as the basic flow. Pursuit of this information will involve three-dimensional coupled 
nonlinear, non-autonomous differential equations. 
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Figure 2. Maximum amplitude of die velocity components as a function of time for Mode I 
disturbances for three values of the z -wavenumber. In each sub-figure, four values of the 
three-dimensionality constant chosen are: a) axisymmetric flow, X = 1 (shown in the figures 
as a small-dashed curve), b) X = 0.5, c) planar stagnation-point flow, X = 0 (shown in the 
figures as a large-dashed curve), and d) X = —0.2. 
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Figure 3. Maximum amplitude of the velocity components as a function of time for Mode II 
disturbances for three values of the z -wavenumber. In each sub-figure, four values of the 
three-dimensionality constant chosen are: a) axisymmelric flow, X = I (shown in the figures 
as a small-dashed curve), b) X = 0.5, c) planar stagnation-point flow, X = 0 (shown in the 
figures as a large-dashed curve), and d) X = -0.2. 
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Figure 4. Normalized energy £ as a function of time. In each sub-figure, four values of the 
three-dimensionality constant chosen arc: a) axisymmetric flow, X = 1 (shown in the figures 
as a small-dashed curve), b) planar stagnation-point fiow, \ = 0, c)\ = -l/3 (shown in the 
figures as a large-dashed curve), and d) X = -0.4. 
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